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ANALYSIS  OF  TURBULENT  FREE-CONVECTION  BOUNDARY  LAYER 

ON  FLAT  PLATE 

By  E.  R.  G.  Eckert  and  Thomas  W.  Jackson 


SUMMARY 

With  the  use  of  Kerman's  integrated  momentum  equation  for  the 
boundary  layer  and  data  on  the  wall -shearing  stress  and  heat  transfer 
in  forced- convection  flow,  a  calculation  was  carried  out  for  the  flow 
and  heat  transfer  in  the  turbulent  free-convection  boundary  layer  on  a 
vertical  flat  plate.  The  calculation  is  for  a  fluid  with  a  Prandtl 
number  that  is  close  to  unity. 

A  formula  was  derived  for  the  heat-transfer  coefficient  that  was 
in  good  agreement  with  experimental  data  in  the  range  of  Grashof  numbers 
from  10  to  101  .  Because  of  the  good  agreement  between  the  theoreti¬ 
cal  formula  and  the  experimental  data  the  formula  may  be  used  to  obtain 
data  for  higher  Grashof  numbers.  The  calculation  also  yielded  formulas 
for  the  maximum  velocity  in  the  boundary  layer  and  for  the  boundary- 
layer  thickness. 


INTRODUCTION 

Recent  developments  in  the  field  of  gas  turbines  have  revealed  the 
need  for  data  on  heat  transfer  in  turbulent  free-convection  flow  at 
very  high  Grashof  numbers  (10^  to  10^).  For  instance,  in  using  the 
method  of  free-convection  cooling  for  turbine  blades,  the  centrifugal 
forces  generate  a  free-convection  flow  of  the  cooling  fluid  in  the 
blade  passages  that  is  within  the  preceding  range  of  Grashof  numbers. 
Free-convection  flow  is  also  superimposed  on  the  flow  of  the  cooling 
air  in  the  hollow  blades  of  air-cooled  gas  turbines  and  may  consider¬ 
ably  influence  the  heat  transfer  under  certain  conditions.  The  radial 
flow  present  in  the  cooled  boundary  layers  oh  the  outside  of  cooled 
turbine  blades  is  also  of  a  type  similar  to  the  free-convection  flow 
in  the  blade  coolant  passages. 

In  order  to  understand  and  evaluate  such  cooling  processes, 
information  on  the  heat  transfer,  on  the  boundary- layer  thickness,  and 
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on  the  velocities  connected  with  free- convection  flow  is  necessary.  The 
knowledge  of  turbulent  free- convection  flow,  however,  is  limited. 
Experiments  on  plane  vertical  surfaces  give  heat-transfer  data  up  to 
Grashof  times  Prandtl  numbers  of  1012  (summarized  in  references  1  and 
2).  Griffiths  and  Davis  (reference  3)  determined  in  addition  tempera¬ 
ture  and  velocity  profiles.  Watzinger  and  Johnson  (reference  4)  meas¬ 
ured  heat  transfer  and  temperature  and  velocity  profiles  in  a  super¬ 
imposed  forced-  and  free- convection  flow  in  a  vertical  tube.  Brown  and 
Marko  (reference  5)  show  by  theoretical  considerations  that  a  general 
relation  exists  between  the  Grashof  number  that  characterizes  free- 
convection  flow  and  the  Beynolds  number  that  characterizes  forced  flow. 
Colburn  and  Hougen  (reference  6)  derive  a  formula  for  the  heat  transfer 
in  turbulent  free-convection  flow  on  a  vertical  plate.  However,  only 
the  laminar  sublayer  in  the  immediate  neighborhood  of  the  wall  was 
investigated  in  reference  6  and  the  thickness  of  this  sublayer,  made 
dimensionless  by  the  shearing  stress  velocity  and  the  kinematic  vis¬ 
cosity,  was  assumed  to  be  the  same  as  in  forced  flow.  This  analysis, 
therefore,  gives  no  information  on  the  whole  boundary- layer  thickness 
and  on  the  velocities. 

The  problem  of  turbulent  free-convection  flow  is  approached  herein 
using  another  method.  Certain  shapes  are  assumed  for  the  temperature 
and  velocity  profiles  in  the  free-convection  boundary  layer.  In  addi¬ 
tion,  an  empirical  relation  for  the  shearing  stress  on  the  wall  and  a 
heat-transfer  coefficient  derived  from  forced- convection  flow  are  used 
to  estimate  the  boundary- layer  thickness,  the  maximum  velocity  within 
the  boundary  layer,  and  a  heat-transfer  coefficient  for  free-convection 
flow.  No  experimental  free-convection  measurements  are  used  in  the 
calculations. 


SHEARING  STRESS  AMD  HEAT  FLOW  IN  FORCED- 
CONVECTION  BOUNDARY  LAYER 

The  relations  for  forced- convection  flow  that  are  needed  in  the 
theoretical  free-convection  calculations  are  compiled  in  this  section. 
It  is  known  that  for  Reynolds  numbers  that  are  not  too  high,  the 
velocity  profile  in  a  turbulent  boundary  layer  on  a  flat  plate  can  be 
represented  by  the  equation5 
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(All  symbols  are  defined  in  the  appendix.)  The  shearing  stress  on  the 
plate  surface  in  such  a  flow  is  given  by  the  equation  (reference  7), 


T 


W 


0.0225  pu2 


(2) 


This  equation  is  the  relation  for  the  shearing  stress  on  the  wall  that 
is  used  in  the  derivation  of  the  boundary- layer  equations  for  the  tur¬ 
bulent  free- convection  boundary  layer.  By  introducing  the  expression 
for  the  velocity  profile,  the  shearing  stress  can  also  be  expressed  by 
the  velocity  u^  at  the  outer  edge  of  the  boundary  layer  and  by  the 

boundary- layer  thickness  5 

3-/4 

Tw  -  0.0225  p%  (Jg)  (3) 


For  fluids  or  gases  with  a  Prandtl  number  equal  to  one,  the  tem¬ 
perature  profile  is  similar  in  shape  to  the  velocity  profile.  When  the 
temperatures  within  the  boundary  layer  are  measured  with  the  temperature 
outside  the  boundary  layer  as  reference,  the  expression  for  the  tem¬ 
perature  profile  is 


*-*v[l-(f)l/7]  <*> 

Equation  (4)  is  valid  for  values  of  y  less  than  6.  For  larger  values 
of  y,  0=0. 

Reynolds  analogy  between  the  turbulent  exchange  of  momentum  and 
heat  gives  the  relation  (reference  8) 


Because  the  temperature  and  velocity  profiles  are  similar  in  shape  for 
a  Prandtl  number  equal  to  one  the  finite  temperature  and  velocity  dif¬ 
ferences  at  two  arbitrary  points  in  the  flow  can  be  introduced  into  the 
last  equation  instead  of  the  differentials.  Using  the  differences 
between  the  wall  and  the  flow  outside  the  boundary  layer  and  specifying 
the  heat  flow  and  the  shearing  stress  for  the  position  at  the  wall 
give  the  following  equation: 


% 

Tv 


6v 


(6) 
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The  introduction  of  the  law  (equation  (3))  for  the  shearing  stress  into 
this  expression  yields  the  following  equation: 


=  0.0225  gpcpu1©w 


(7) 


Heat  transfer  is  usually  calculated  with  dimensionless  moduli.  Such  a 
modulus  is  the  Stanton  number 


St 


% 

gPcpul®v 


Experimental  investigations  show  that  for  fluids  with  Prandtl  numbers 
from  approximately  0.5  to  50  (reference  1,  p.  520)  the  same  relation 
holds  when  it  is  multiplied  by  the  factor  Pr“2/3. 


St 


% 

gpCpUi^w 


(8) 


This  relation  for  the  heat  flow  qv  is  used  in  the  derivation  of  the 
boundary- layer  equations  for  the  turbulent  free- convection  boundary 
layer. 


Replacing  the  boundary- layer  thickness  5  by 


5 


which  is  the  boundary- layer  thickness  on  a  flat  plate  in  turbulent 
forced- convection  flow  (reference  1,  p.  481),  transforms  equation  (8) 
into  the  widely  accepted  formula  for  heat  transfer  on  a  flat  plate  in 
the  turbulent  range 


Nu 


(9) 


DETERMINATION  OF  TURBULENT  FREE- CONVECTION  BOUNDARY  LAYER 
Derivation  of  Boundary- Layer  Equations 

k 

If  a  stationary  plane  vertical  wall  is  heated  to  a  temperature 
higher  than  the  surroundings,  the  layer  of  fluid  adjacent  to  the  wall 
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is  heated  by  conduction  from  the  wall.  In  this  way,  buoyancy  forces 
are  generated  that  cause  this  layer  to  flow  in  an  upward  direction. 

This  layer  of  fluid  adjacent  to  the  wall  to  which  the  vertical  motion 
is  confined  is-  called  the  free- convection  boundary  layer.  The  boundary 
layer  begins  with  zero  thickness  at  the  lower  end  of  the  vertical  wall 
and  increases  in  thickness  in  the  upward  direction.  The  boundary  layer 
becomes  turbulent,  depending  on  the  critical  Grashof  number,  at  a  cer¬ 
tain  distance  from  the  lower  end  of  the  wall.  The  distance  measured 
vertically  from  the  lower  end  of  the  wall  is  called  x  and  the  dis¬ 
tance  normal  to  the  wall  y.  In  order  to  determine  the  boundary- 
layer  thickness  for  steady  state,  a  small  stationary  volume  element 
in  the  turbulent  region  of  the  boundary  layer  is  considered.  Figure  1 
shows  this  volume  element.  The  dimensions  of  the  element  are  dx  along 
the  wall  and  Z  normal  to  the  wall.  The  length  Z  should  be  larger 
than  the  boundary-layer  thickness  6.  For  two-dimensional  flow,  the 
dimension  of  the  volume  element  normal  to  the  plane  of  figure  1  may  be 
considered  to  be  unity.  The  upward  velocity  of  the  fluid  in  plane  1-1 
at  a  distance  y  from  the  surface  of  the  wall  is  u.  Then  the  mass 
flow  through  a  small  area  with  a  width  dy  is  pu  dy  and  the  flow  of 
momentum  in  the  x-direction  is  pu2dy.  The  momentum  flow  in  the 
x-direction  entering  the  volume  element  through  plane  1-1  is 

J\ u2dy 

In  progressing  to  plane  2-2  the  momentum  flow  changes  by 


The  mass  flow  entering  the  volume  element  through  plane  1-1  is 
generally  different  from  the  mass  flow  leaving  the  element  through 
plane  2-2.  Therefore  fluid  enters  or  leaves  the  volume  element  through 
the  plane  parallel  to  the  wall  at  a  distance  Z .  Because  it  is  assumed 
that  the  velocity  in  the  x-direction  outside  the  boundary  layer  is  so 
small  that  it  can  be  neglected,  no  momentum  in  the  x-direction  is 
carried  through  the  plane. 

.  The  rate  of  change  of  momentum  must  be  in  equilibrium  with  the 
forces  acting  in  the  x-direction  on  the  fluid  within  or  on  the  surface 
of  the  volume  element  considered.  A  shearing  stress  Ty  acts  on  the 

wall.  The  force  connected  with  this  stress  is  Tv  dx.  No  shearing 
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stress  occurs  on  the  surface  of  the  volume  element  that  is  parallel  to 
the  vail  and  at  a  distance  l  from  the  wall  because  outside  the  bound¬ 
ary  layer  the  velocity  in  the  x- direction  is  zero. 

According  to  the  boundary- layer  theory,  the  pressure  change  can  be 
neglected  along  any  normal  to  the  surface.  A  constant  pressure  differ¬ 
ence  dp  therefore  exists  between  planes  1-1  and  2-2.  This  pressure 
difference  gives  a  force  on  the  volume  element  of  magnitude  2  dp.  In 
addition  to  the  forces  on  the  surfaces  there  is  a  force  due  to  the 
weight  of  the  fluid  within  the  volume  element 


The  process  of  summing  up  all  the  forces  and  equating  them  to  the 
change  in  momentum  flow  gives  the  momentum  equation 


No  flow  exists  outside  the  boundary  layer  and  therefore  the  pres¬ 
sure  difference  between  planes  1-1  and  2-2  is  balanced  by  the  weight  of 
the  fluid  layer  between  the  planes 


dp  =  p6g  dx 


When  both  sides  of  this  equation  are  multiplied  by  2  and  the  right- 
hand  side  of  the  equation  is  changed  to  the  integral  form,  the  follow¬ 
ing  equation  is  obtained 


Id p  = 


Introducing  the  preceding  equation  into  the  momentum  equation  gives 


pu^  dy  dx  =  g 


(p5-p)  d7 


0 


dx  -  T  dx 
w 


A 
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Introducing  the  expansion  coefficient  defined  by  the  equation 


the  first  term  on  the  right-hand  side  of  the  preceding  equation  can  be 
transformed  into 


Pp(t-te)  dy 


Designating  the  difference  between  the  temperature 


and  the  temperature  t§ 
expression  is  obtained: 


outside  the  boundary  layer  by 


at  the  distance  y 
6 ,  the  following 


g  dx  pp0  dy 

Jo 

In  the  applications  considered,  the  essential  influence  of  the 
density  changes  on  the  flow  is  taken  into  account  by  the  introduction 
of  the  expansion  coefficient  (3.  The  other  influences  of  variable  density 
on  the  flow  and  the  variation  of  the  expansion  coefficient  ,6  with 
small  temperature  differences  are  negligible.  Therefore,  p  and  p 
can  be  assumed  constant  in  the  preceding  expressions.  The  momentum 
equation  for  the  free- convection  boundary  layer  therefore 
becomes 


u2  dy 


0  dy  -T 


A  similar  equation  is  set  up  for  the  heat  flow  through  the  volume 
element  in  figure  1.  The  heat  carried  with  the  fluid  through  plane  1-1 


u  6  dy 


where  the  enthalpies  Cp0  are  measured  from  the  temperature  outside 
the  boundary  layer.  The  specific  heat  and  density  are  considered 
constant . 
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The  heat  carried  out  of  plane  2- 
plane  1-1  by 


gpc 


2  differs  from  that  carried  into 


u0  dy  /  dx 


The  difference  in  the  heat  flow  through  planes  1-1  and  2-2  must  come 
from  the  surface  and  the  heat  flow  leaving  the  plate  per  unit  time  and 
area  is  therefore  y 


^w 


(id 


Equation  (ll)  is  the  heat-flow  equation  for  the  free- convection  boundary 
layer . 


Equations  for  the  laminar  free- convection  boundary  layer  that  are 
analogous  to  equations  (10)  and.  (ll)  were  derived  and  solved  by  Squire. 
(See  reference  9.) 


Solution  of  Boundary-Layer  Equations 

Equations  (10)  and  (ll)  are  sufficient  to  calculate  the  boundary- 
layer  thickness  6  and  the  velocity  u  when  the  shape  of  the  velocity 
and  temperature  profiles  within  the  boundary  layer  and  the  laws  for  the 
shearing  stress  and  the  heat  flow  on  the  wall  are  known.  It  was  sug¬ 
gested  by  von  Karman  that  an  approximate  solution  may  be  obtained  by 
introducing  approximate  shapes  for  both  profiles  (reference  7).  The 
accuracy  of  the  solution  is  better  the  more  closely  the  assumed  shapes 
correspond  to  the  real  profiles.  The  method,  however,  proved  compar¬ 
atively  insensitive  to  changes  in  the  profile  shape. 

As  mentioned  in  the  section  " INTRODUCTION"  some  information  on  the 
shape  of  the  velocity  and  temperature  profiles  in  turbulent  free- 
convection  flow  can  be  obtained  from  reference  3 .  In  figure  2 ,  the 
measured  values  of  both  profiles  in  the  turbulent  range  are  presented 
with  the  ratios  ©/©w  and  u/^  plotted  against  the  distance  from 

the  wall  in  an  arbitrary  scale.  The  value  uDiax  was  determined  by 

drawing  curves  through  the  measured  velocity  values  given  in  table  VI 
of  reference  3. 
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The  distance  from  the  wall  where  the  curves  u/u^  have  the  value 
0.5  was  arbitrarily  assigned  the  value  1.  The  lines  in  figure  2  are  the 
curves  by  which  the  temperature  and  velocity  profiles  are  approximated 
herein.  These  profiles  are  presented  in  figure  2  in  such  a  way  that  the 
velocity  profile  has  the  value  0.5  at  a  distance  1  and  the  temperature 
profile  the  value  0.2  at  a  distance  0.5  from  the  wall.  An  indication  of 
the  shape  of  turbulent  free- convection  flow  profiles  can  also  be  obtained 
from  the  measurements  by  Watzinger  and  Johnson  (reference  4)  on  mixed 
forced-  and  free- convection  flow  in  a  vertical  tube.  Some  of  these 
measured  profiles  are  presented  in  figure  3. 

Equation  (4),  which  describes  the .forced- convection  temperature 
profile,  is  used  for  the  temperature  profile  in  turbulent  free- convection 
flow.  The  profiles  of  figure  2  indicate  that  this  assumption  is  rea¬ 
sonable.  The  velocity  profile  in  free- convection  flow,  however,  differs 
from  the  forced-flow  profile  by  having  a  velocity  of  zero  outside  the 
boundary  layer.  Accordingly,  equation  (l)  is  multiplied  by  a  function 
of  y/6,  which  brings  the  velocities  back  to  the  value  zero  at  the 
outer  edge  of  the  boundary  layer  where  y  =  6.  Two  simple  expressions 
that  fulfill  this  condition  and  that  agree  quite  well  with  the  meas¬ 
ured  values  are  plotted  in  figure  2.  The  velocity  profile  that  rep¬ 
resents  the  measured  points  somewhat  better  is  used  in  the  following 
calculations : 


(12) 


Equation  (12 )  is  valid  for  values  of  y  less  than  6.  For  larger 
values  of  y,  u  =  0.  The  maximum  velocity  umay  of  this  profile  can 
be  found  by  differentiating  equation  (12),  giving 


u. 


max 


=  0.537  u-j_ 


(13) 


By  use  of  equations  (4)  and  (12)  and  the  fact  that  u  =  0  and  0=0 
for  distances  greater  than  6  from  the  wall,  the  integrals  in  the 
momentum  and  heat-flow  equations  (equations  (10)  and  (ll))  become 


I 

/  u2  dy  =  0.0523  bu±2 

J  0 


0.125 

w 
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The  same  value  of  5  was  used  for  the  velocity  and  temperature 
profiles.  Calculations  on  forced- convection  heat  transfer  indicate 
that  the  same  value  of  6  for  the  velocity  and  temperature  profiles 
can  be  used  for  fluids  that  have  a  Prandtl  number  close  to  one  (ref¬ 
erence  8).  For  very  large  or  small  Prandtl  numbers  in  forced  flow, 
the  thickness  of  the  temperature  boundary  layer  is  considerably  smaller 
or  larger  than  the  thickness  of  the  velocity  boundary  layer.  This  con¬ 
dition  is  also  probably  true  for  free- convection  flow. 

In  order  to  solve  boundary- layer  equations  (10)  and  (11)  the  laws 
for  the  shearing  stress  on  the  wall  and  the  heat  flow  must  be  intro¬ 
duced.  It  is  assumed  that,  in  the  layers  very  near  the  wall,  the 
conditions  are  similar  for  free-  and  forced-convection  flow  and  that 
the  same  laws  for  the  shearing  stress  and  the  heat  flow  that  are  used 
in  forced- convection  flow  can  therefore  be  used  in  free- convection 
flow.  A  similar  assumption  for  the  shearing  stress  was  used  by  von 
Kerman  to  calculate  the  radial  flow  on  a  rotating  disk,  with  satis¬ 
factory  results  (reference  7). 

As  a  consequence  of  the  preceding  assumption  equations  (2)  and 
(8)  will  be  used  for  free-convection  flow.  Very  near  the  wall  the 
second  term  in  equation  (12)  for  the  free-convection-velocity  profile 
is  one  and,  therefore,  the  same  transformation  that  led  in  forced  flow 
from  equation  (2)  to  equation  (3)  can  be  made  for  free-convection  flow. 
By  introducing  the  evaluated  integrals  and  equations  (3)  and  (8)  the 
momentum  and  heat-flow  equations  become 


0.0523  ~- 
dx 


0.125  gpe^s 


0.0225  u. 


0.0366 


(v)  ■ 


0.0225  u. 


Equations  (14)  and  (15)  are  total  differential  equations  from  which  the 
two  unknown  values  5  and  U]_  may  he  determined  as  functions  of  x. 

Equations  (14)  and  (15),  can  be  solved  by  introducing 


ui  =  cux 


(16) 
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8  =  C=x 


n 


(17) 


The  introduction  of  these  values  transforms  equations  (14)  and  (15)  into 


0.0523 (2m+n)  Cu2C&x2m+n-1  =  O.125gO0wC8xn 


0.0225  CU7A  C6‘1''4  V1/4  x(7»A)-(n/4) 


(18) 


0.0366(m+n)CuC5x 


m+n-1 


=  0.0225  (Pr)"  2/3  vl/4  Cu3/4  C8-l/4x(3»A)-(n/4)  (19) 

Because  equations  (18)  and  (19)  must  be  valid  for  any  value  of  x,  the 
exponents  of  x  must  be  identical,  that  is, 


2m+n-l  =  n  • 

m+n-1  = 


7m  n 
4  “  4 

3m  n 


It  can  be  seen  that  this  condition  is  fulfilled  if 

m  =  i  and  n  =  ^ 


(20) 


With  these  values  of  m  and  n  equation  (19)  can  be  solved  for  the 
constant  Cu 

Cu  =  0.0689  VC5"5(Pr)“8/3 
Introducing  this  constant  into  equation  (18)  gives 

CB10  =  0.00338  -|g-  [l+0.494(Pr)2/3]  (Pr)"16/3 

gpe^x0 

Introducing  the  Grashof  number  Gr  =  -  gives  for  the  velocity 


(21) 


(22) 


ux  =  1.185  *  (Gr)1//2 £l+0. 494(Pr)2//3] 


-1/2 


(23) 
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and  for  the  boundary-layer  thickness 


6  =  0.565  x(Gr)"1/10(pr)-8/15  |jL+0.494(Pr)2/3] 


1/10 


(24) 


More  significant  than  the  value  is  the  maximum  velocity 

umax  boundary  layer.  (See  equation  (13).)  A  Reynolds 

number  is  obtained  to  make  this  velocity  dimensionless 


-1/2 


Ee 


max 


=  =  0.636  (Gr)1/2  [jL+0.494(Pr)2/3 J  (25) 


A  more  significant  value  characterizing  the  thickness  of  the  bound¬ 
ary  layer  than  the  value  6  used  up  to  now  is  the  displacement  thick¬ 
ness  5*,  which  is  defined  for  free-convection  flow  by  the  equation 


6* 


p*  00 


u 

0  %iax 


dy 


With  the  use  of  equations  (12)  and  (13)  it  is  found  that 

6*  =  0.272  5 


The  displacement  thickness  made  dimensionless  by  the  distance  x 
from  the  leading  edge  of  the  plate  is  therefore 

^  =  0.154  (Gr)-1A°(pr)'’8/15  [\+0.494(Pr)2/3J 

X 

The  heat- transfer  coefficient  H  =  qw/$v  at  the  point  x  on  the 

wall  is  found  by  introducing  the  boundary- layer  thickness  6  and  the 
velocity  u1  into  equation  (8).  Changing  to  the  usually  presented 
Nusselt  number  Nu  =  Hx/k  gives 


Nu  =  0.0295 


(Gr)2/5  (Pr)7/15 


1+0.494  (Pr) 


2/3 


n-2/5 


(26) 


In  order  to  compare  this  equation  with  experimental  results,  it  is  nec¬ 
essary  to  change  from  the  local  heat-transfer  coefficient  to  the  average 
value  along  the  plate.  By  introducing  the  expression  for  the  Grashof 
number  into  equation  (26),  it  can  be  seen  that  the  local  heat- transfer 
coefficient  is  proportional  to  the  power  0.2  of  the  distance  x 
(H  =  CgX0,2).  With  the  assumption  that  the  boundary  layer  is  turbulent 
from  the  leading  edge,  the  average  heat- transfer  coefficient  becomes 
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In  reality  the  boundary  layer  is  first  laminar  and  only  at  a  certain 
distance  from  the  lower  edge  of  the  plate  does  it  become  turbulent. 

The  preceding  expression  for  the  average  heat- transfer  coefficient  can 
therefore  be  expected  to  represent  the  true  value  only  at  Grashof  num¬ 
bers  so  high  that  the  extent  of  the  laminar  boundary  layer  is  small 
compared  with  the  total  length  x.  This  limit  for  the  Grashof  number 
seems  to  be  near  10^®.  For  Grashof  numbers  higher  than  this  value  the 
average  Nusselt  number  can  be  calculated  with  equation  (27);  therefore. 


Nu 


=  0.0246  (Gr)2/5  (Pr)7/15  jl+0.494  (Pr)2/3 


av 


-2/5 


(28) 


In  order  to  determine  to  what  extent  the  shape  of  the  velocity 
profile  influences  the  results,  a  calculation  was  made  with  the  second 
velocity  distribution  shown  in  figure  2* 


(29) 


Only  the  numerical  constants  are  influenced  by  the  change  in  the  profile 
shape.  The  constant  preceding  the  value  (Pr)2/3  changes  from  0.494  to 
0.342  in  all  equations.  In  addition,  the  value  0.636  changes  to  0.487 
in  the  equation  for  Remax;  in  the  equation  for  S*/x  the  value  0.154 
changes  to  0.200  and  in  the  equation  for  Nuav  the  constant  0.0246 
changes  to  0.0198. 


COMPARISON  WITH  EXPERIMENTS 

In  figure  4  the  results  of  experiments  carried  out  in  different 
investigations  (references  1  and  2)  are  plotted  as  the  average  Nusselt 
number  Nuav  against  the  product  Gr  Pr.  For  the  lower  values  of 
Gr  Pr  the  experimental  results  quite  accurately  fit  the  equation 

Nuav  =  0.555  (Gr  Pr)1/4 

which  was,  with  a  slight  adjustment  of  the  constant,  theoretically 
derived  for  laminar  free- convection  flow.  The  experimental  results  in 
the  turbulent  range  (Gr  Pr  =  1010  to  1012)  may  be  represented  by 
the  equation 
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Wuav  =  0.0210  (Gr  Pr)2/5  (30) 

The  exponent  of  the  Grashof  number  in  this  equation  is  the  same  as  that 
of  equation  (28)  derived  in  the  previous  section.  When  equation  (28) 
is  transformed  into  the  form  of  equation  (30 )  in  such  a  way  that  the 
values  of  both  equations  are  the  same  for  Pr  =  0.72,  the  constant  of 
equation  (28)  becomes  0.0210.  The  heat- transfer  coefficient  derived 
with  equation  (12)  for  the  velocity  profile  is  in  perfect  agreement 
with  the  experimental  results.  Such  agreement  is  probably  a  coinci¬ 
dence.  The  profile  given  by  equation  (29),  which  does  not  fit  the 
velocity  distributions  measured  by  Griffiths  and  Davis  in  reference  3 
as  well  as  does  the  first  profile,  gives  heat-transfer  coefficients 
that  are  17-percent  lower  than  the  measured  values. 

The  values  for  the  maximum  velocity  within  the  boundary  layer  and 
the  boundary-layer- displacement  thickness  agree  poorly  with  the  values 
measured  in  reference  3  in* the  turbulent  range.  Heat-transf er  coeffi¬ 
cients  were  not  measured  therein  on  the  experimental  apparatus  on  which 
velocity  and  temperature  profiles  were  obtained. 

Whereas  the  velocity  and  temperature  profile  shapes  as  measured  in 
reference  3  are  typical  of  turbulent  free-convection  flow  the  order  of 
magnitude  of  these  profiles  appears  to  be  in  error.  It  can  be  shown 
that  there  is  disagreement  within  the  measured  values  themselves.  The 
heat  given  off  by  the  plate  to  the  air  stream  must  be  carried  away 
within  the  boundary  layer.  The  measured  temperature  and  velocity  pro¬ 
files  as  well  as  the  measured  maximum  velocity  and  the  boundary- layer 
thickness  vary  little  along  the  plate  in  the  turbulent  range ,  which 
means  that  only  a  small  part  of  the  heat  given  off  by  the  wall  is  found 
in  the  boundary  layer.  The  horizontal  dimension  of  the  plate  may  have 
been  too  small  compared  with  the  vertical  length  to  make  the  flow  in 
the  center  part  two-dimensional  and  air  may  have  flowed  into  the  bound¬ 
ary  layer  from  the  sides. 


SUMMABY  OF  RESULTS 

With  the  use  of  Kerman1  s  approximate  method,  a  calculation  was 
carried  out  for  the  flow  and  heat  transfer  in  the  turbulent  free- 
convection  boundary  layer  on  a  vertical  flat  plate.  The  calculation 
used  relations  for  the  heat  flow  and  shearing  stress  on  the  wall  devel¬ 
oped  for  forced  turbulent  flow  and  velocity  and  temperature  profiles 
that  approximate  well  the  shapes  of  profiles  measured  by  Griffiths  and 
Davis. 
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A  formula  was  derived  for  the  heat- transfer  coefficient  that  was  in 
good  agreement  with  measured  values  in  the  range  of  Grashof  numbers 
from  lO-'-O  to  10^  and  that  can  be  used  to  extrapolate  the  values  into 
the  range  of  higher  Grashof  numbers.  The  formula  is  valid  for  Prandtl 
numbers  that  are  close  to  unity. 

The  calculation  also  yielded  formulas  for  the  maximum  velocity  in 
the  boundary  layer  and  for  the  boundary- layer  thickness. 


Lewis  Flight  Propulsion  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  July  12,  1950. 
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APPENDIX  -  SYMBOLS 

The  following  symbols  are  used  in  this  report: 
thermal  diffusivity,  sq  ft/sec 

0  9 

constant  for  variation  of  heat- transfer  coefficient  with  x 

constant  for  variation  of  velocity  u  in  boundary  layer 
with  xm 

constant  for  variation  of  boundary- layer  thickness  with  xn 
specific  heat  at  constant  pressure,  Btu/(lb)(°F) 


Grashof  number. 


V 


acceleration  due  to  gravity,  ft/sec^ 
heat-transfer  coefficient,  Btu/(sq  ft)(sec)(°F) 
heat  conductivity,  Btu/(ft)(sec)(°F) 
length  (fig.  1),  ft 
exponent 

Nusselt  number,  — 


Pr 

P 

q 

■Remax 

St 


exponent 

Prandtl  number, 

'  a  k 

pressure,  lb/(sq  ft) 

specific  heat  flow,  Btu/(sec)(sq  ft) 

Reynolds  number  based  on  maximum  velocity 

1/4 


u  x 
max 

V 


Stanton  number. 


*w 


gPGpUl®W 


=  0.0225 


&) 
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t  temperature,  °F 

u  velocity  component  in  x- direct ion,  ft/sec 

u*^  velocity  outside  boundary  layer  of  comparable  forced- convection 

flow,  ft/sec 

v  specific  volume,  cu  ft/lb 

x  coordinate  (distance  along  plate  from  starting  point  of  bound¬ 

ary  layer),  ft 

y  coordinate  (distance  from  wall),  ft 

(3  expansion  coefficient,  l/°F 

6  boundary-layer  thickness,  ft 

6*  displacement  thickness  of  boundary  layer,  ft 

6  temperature  difference,  °F 

Qv  temperature  difference  between  wall  and  fluid  outside  of  bound¬ 

ary  layer,  °F 

M-  absolute  viscosity,  pv>  lb  sec/sq  ft 

V  kinematic  viscosity,  sq  ft /sec 

p  mass  density,  (lb)(sec2)/ft^ 

T  shearing  stress,  lb/sq  ft 

Subscripts: 
av  average  value 

max  maximum  value 

w  on  wall 
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Figure  1*  -  Nomenclature  in  deriving  equations  for  free -convection  "boundary  layer. 
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Figure  2.  -  Experimental  and  theoretical  velocity  and  temperature  profiles  for  turbulent 
free-convection  flow  on  a  flat  plate.  Experimental  data  from  reference  3. 
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Figure  4.  -  Average  Nusselt  number  for  free- convection  flows  on  a  vertical  plate 
Experimental  data  from  references  1  and  2. 


